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ABSTRACT 

CEINCI-LAB is a computer system that allows static or dynamic analysis of 
structures. This article presents the interaction of CEINCI-LAB functions with 
OpenSees commands for the development of nonlinear static analysis using the 
Pushover technique of a steel structure. The steel structure considered in this 
study is a 6-story, with 4 column axes, and different structural reinforcement 
mechanisms such as the chevron brace and multi-story X-brace. 

OpenSees requires a certain effort to enter the geometry of the frames of 
the structures, involving more time when working multi-story structures, since each 
of the nodes that will make up the elements must be defined, be they beams, 
columns or diagonals in the case of structures with reinforcement and indicate how 
with these nodes each element is formed with a particular property and material. 
This article seeks to emphasize the benefits of developing analysis with the help of 
these two computational tools. The new functions of CEINCI-LAB facilitate the 
entry of multi-story structures and the creation of nodes in the center of the 
openings to analyze the reinforced structure. 

Both CEINCI-LAB and OpenSees can generate this analysis (Pushover) to 
form directly without the interaction between them, however, by coupling certain 
CEINCI-LAB functions, data entry and presentation of results is facilitated and the 
user can focus on understand the behavior of the structure without the complexity 
that programming in OpenSees can represent. 

keywords: OpenSees. Computer system CEINCI-LAB. Nonlinear Static Analysis. 

Pushover. Seismic Reinforcement with steel braces. 
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ACOPLAMIENTO DE OPENSEES CON CEINCI-LAB 
PARA ANALISIS ESTATICO NO LINEAL. 
PRIMERA PARTE: REFORZAMIENTO SISMICO CON 
DIAGONALES DE ACERO 

RESUMEN 

CEINCI-LAB es un sistema de computacion que permite el analisis 
estatico o dinamico de estructuras. En este articulo se presenta la interaccion de 
las funciones de CEINCI-LAB con comandos de OpenSees para el desarrollo del 
analisis estatico no lineal mediante la tecnica de Pushover de una estructura de 
acero. La estructura en analisis es de 6 pisos, con 4 ejes de columnas, ademas se 
presenta diferentes mecanismos de reforzamiento estructural en base a 
diagonales de acero ubicadas en forma de V invertidas y de X de piso a piso. 

OpenSees requiere cierto esfuerzo para el ingreso de la geometrfa de los 
porticos de las estructuras, involucrando mayor tiempo cuando trabajamos 
estructuras de varios pisos, ya que se debe definir cada uno de los nudos que 
conformaran los elementos, sean estos vigas, columnas o diagonales en el caso 
de estructuras con reforzamiento e indicar como con estos nudos se forma cada 
elemento con una propiedad y material particular. Este artfculo busca enfatizar las 
bondades de desarrollar el analisis con la ayuda de estas dos herramientas 
computacionales. Las nuevas funciones de CEINCI-LAB facilitan el ingreso de 
estructuras de varios pisos y la creacion de nudos en el centro de los vanos para 
analizar el portico reforzado. 

Tanto CEINCI-LAB como OpenSees pueden generar este analisis 
(Pushover) de formar directa sin la interaccion entre ellos, sin embargo, al acoplar 
ciertas funciones de CEINCI-LAB el ingreso de datos y la presentacion de 
resultados se facilita y el usuario puede enfocarse en entender el comportamiento 
de la estructura sin la complejidad que puede representar la programacion en 
OpenSees. 

Palabras Claves: OpenSees. Sistema de computacion CEINCI-LAB. Analisis 
Estatico No Lineal. Pushover. Reforzamiento con diagonales de acero. 

1. INTRODUCCION 

La Tecnica del Pushover o Tecnica del Empujon, tuvo un gran desarrollo a 
finales del siglo XX y todavfa se continua trabajando en esta tematica, porque 
ayuda a visualizar el probable desempeno que va a tener una estructura durante 
un terremoto. Importantes investigaciones se han realizado aplicando programas 
como DRAIN-2DX, para analizar conexiones en estructuras de hormigon armado 
utilizando analisis estatico no lineal, que es otra forma de conocer al Pushover. 
Godfnez y Lopez Almansa, (2014). 

Se ha estudiado la disipacion de estructuras de acero con diagonales 
concentricas y excentricas; mediante el analisis dinamico incremental que se 
fundamenta en el Pushover. E. Del Valle et al. (2014). De igual manera el efecto 
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de los modos superiores en estructuras de hormigon con contraventeo metalico 
tipo Chevron, mediante un Pushover modal. Godinez y Tena, (2014). 

Experimentalmente se ha validado la curva de capacidad sismica que 
relaciona el cortante basal con el desplazamiento lateral maximo, encontrado con 
la tecnica del Pushover y se vio que era muy similar a la que se obtuvo al ensayar 
varios porticos de hormigon armado de un vano y un piso ante cargas laterales 
monotonicas. Diaz et a/. (2015). Se ha presentado unicamente 3 referencias sobre 
las aplicaciones que se tienen a partir de los resultados que se hallan en base al 
Pushover, y una validacion experimental. Con esto lo que se desea mostrar es que 
este tema es muy actual y es asf como se han realizado valiosos aportes al 
sistema de computacion CEINCI-LAB para hallar la curva de capacidad slsmica 
resistente. Para citar un caso se hace referencia al trabajo de Mora y Aguiar (2015, 
2) donde se presenta la Tecnica del Pushover Multimodal, aplicado al analisis de 
estructuras con disipadores de energfa ADAS o TADAS. 

Por otra parte, se han hecho importantes contribuciones a la Ingenierfa 
Sfsmica, por medio de la Tecnica del Pushover en el analisis de estructuras de 
acero, utilizando el sistema OpenSees, Open System for Earthquake Engineering 
Simulation (Tapia y Gutierrez, 2016). 

En la figura 1a, se tiene una estructura de acero de 6 pisos, la misma que 
puede ser reforzada con diagonales en forma de “V” invertida (figura 1b), donde en 
cada piso una diagonal trabaja a compresion y la otra a traccion; se debe tener 
cuidado que la diagonal a compresion no ingrese al rango no lineal, ya que se 
puede presentar un comportamiento indeseable en la viga a la que conectan, 
aunque esto se puede evitar utilizando secciones grandes en las diagonales para 
que el factor de esbeltez sea bajo, con lo que se minimiza el problema de pandeo. 

Otra opcion en la que se elimina el problema de pandeo, es reforzar con 
diagonales en forma de “X”, ahf las diagonales en un sentido (piso inferior y piso 
superior) trabajan a compresion y las diagonales en el otro sentido (piso inferior y 
piso superior) trabajan a tension, de tal manera que existe equilibrio en el nudo 
central de la viga; cuando estos elementos desarrollen su maxima capacida a 
compresion y a traccion, se anula y no hay el problema que puede darse en el 
caso anterior. Lo que interesa ilustrar por ahora en la figura 2, es la sencillez con 
que se ingresan los datos al sistema CEINCI-LAB para dibujar las estructuras. 

En la figura 1b, se presenta la estructura reforzada con diagonales en 
forma de “V” invertida y en la figura 1c, se tiene la estructura reforzada con 
diagonales en forma de “X”, de piso a piso. El programa que dibuja las tres 
estructuras se indica a continuacion y se puede ver que es extremadamente 
sencillo dibujarlas con CEINCI-LAB. 

Con los datos sv, sp y con las funciones: geometria_nudo_viga; 
glinea_portico2; gn_portico2, se dibuja la estructura sin diagonales (ver figura 
1a), lo importante es que ya estan definidas las coordenadas X, Y de la estructura 
y los vectores Nl, NJ, que contienen a los nudos iniciales y finales de los 
elementos. 
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Figura 
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a) Estructura sin reforzamiento; b) Estructura reforzada con diagonales 
en forma de “V” invertida; c) Estructura reforzada con diagonales en 
forma de “X”. 
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Codigo para generar la geometri'a de las 
estructuras: 


% E j g eneracion porticos 

%% Primer Portico es una estructura de acero sin reforzamiento 

sv=[7.315;7.315;7.315]; % Luces de los vanos 

sp=[5.486;3.810;3.810;3.810;3.810;3.810]; % Alturas de cada piso 

% Rutinas para generar la geometrfa del portico sin diagonales 
[nv,np,nudt,nudcol,nudvg,nod,nr]=geometria_nudo_viga(sv,sp); 
[X,Y]=glinea_portico2(nv,np,sv,sp,nod,nr); 

[NI,NJ]=gn_portico2(nr, nv, nudt, nudcol, nudvg); 
figure (1) 

dibujoplano(X,Y,NI,NJ) 

%% Segundo Portico se refuerza el vano central con diagonales en forma de 
V invertida 

mar=[2;5;8;11 ;14;17]; % Numero de vano con reforzamiento 

% Rutinas para ingresar diagonales V invertida 
[GEN]=geometria_nudo_diagonales(nv, np, nudt, mar); 
[NI2,NJ2]=gn_portico(GEN); 

[NI,NJ]=gn_portico3(NI,NJ, NI2,NJ2); 
figure (2) 

dibujoplano(X,Y,NI,NJ) 

%% Tercer Portico se refuerza el vano central con diagonales en forma de "X" 

[nv,np,nudt,nudcol,nudvg,nod,nr]=geometria_nudo_viga(sv,sp); 

[X,Y]=glinea_portico2(nv,np,sv,sp,nod,nr); 

[NI,NJ]=gn_portico2(nr, nv, nudt, nudcol, nudvg); 
mar=[2;8;14]; % Numero de vano con "V" invertida 
marv=[5;11 ;17]; % Numero de vano con "V" 

% Rutinas para ingresar diagonales en forma de X 
[GEN]=geometria2_nudo_diagonales(nv, np, nudt, mar, marv); 
[NI2,NJ2]=gn_portico(GEN); 

[NI,NJ]=gn_portico3(NI,NJ, NI2,NJ2); 
figure (3) 

dibujoplano(X,Y,NI,NJ) 


Figura 2 Codigo para generar las geometrfas de las estructuras en analisis. 

Para dibujar las diagonales se debe ver la numeracion de los vanos que se 
tiene en la figura 3, y unicamente en el vector mar se identifican los vanos donde 
se encuentran las diagonales en forma de “V” invertida; lo que resta ahora es 
identificar el nudo inicial y final de cada diagonal, esto se lo hace con las 
funciones: geometria_nudo_diagonales y gn_portico. Finalmente, se incorporan 
los nudos iniciales de las diagonales NI2, NJ2 a los nudos iniciales y finales de la 
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estructura sin diagonales Nl, NJ; la incorporacion se la hace con la funcion 

gn_portico3 

Un procedimiento similar se realiza para colocar las diagonales en forma 
de “X”, primero en el vector mar se indica los vanos de las diagonales en forma de 
“V” invertida y con el vector marv se identifican los vanos con diagonales en “V”. 
La obtencion de los vectores con nudos iniciales y finales de las diagonales en 
forma de “X” se realiza con las funciones: geometria2_nudo_diagonales y con la 
funcion gn_portico. Por ultimo, la incorporacion de estos vectores, de las 
diagonales a los vectores de nudo inicial y final de los elementos sin diagonales se 
lo hace con la funcion gn_portico3. 
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Figura 3 Numeracion de vanos para ingreso de diagonales. 

Con 4 vectores de datos: sv, sp, mar, marv, son suficientes para generar 
los 3 porticos de las figuras 1 a, b, c. Claro esta que esto se realiza con la 
utilizacion de las funciones de CEINCI-LAB, pero realmente son muy pocas. 
(Aguiar et al. 2019) 

Tanto CEINCI-LAB como OpenSees pueden generar este analisis 
(Pushover) de formar directa sin la interaccion entre ellos, sin embargo, por 
facilidad para el ingreso de datos y la presentacion de resultados se utilizan 
funciones de CEINCI-LAB y debido al potencial de comandos desarrollados en 
OpenSees se emplean funciones de este que permiten simular estructuras 
considerando la no linealidad de materiales con modelos que se asemejan de 
mejor manera a la realidad. 

En lo que sigue de este artfculo se va a encontrar la curva de capacidad 
sismica resistente de una estructura de acero y su reforzamiento, utilizando las 
funciones de CEINCI-LAB para la entrada de datos y salida de resultados, 
ademas del calculo propiamente dicho del Pushover que se lo hace con 
OpenSees. 
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2. DESCRIPCION DE LA ESTRUCTURA SIN REFUERZO 

Lo principal es definir la geometrfa del portico en analisis, para este caso 
se tiene una estructura de acero con 3 vanos de 7.315 m y 6 pisos, la altura del 
primer nivel es de 5.480 m mientras que en los niveles superiores son de 3.810 m 
de altura (Fig. 4). 
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Figura 4 Geometrfa de la estructura de acero de 6 pisos e identificacion de las 
secciones. 

Se utilizaron las siguientes propiedades de material para modelar el 
portico: ASTM A36, Fy = 36 ksi, Fu = 40 ksi y E = 29000 Ksi. Las dimensiones de 
las columnas y vigas se definen en la Tabla 1, donde la identificacion del color es 
consistente con la figura 4. 

Tabla 1 Descripcion de las secciones de los elementos estructurales. 


i 



Color Elemento 


Columnas 

W30x173 


Vigas 

W30x99 


Dimensiones (mm) 
d x b x tf x tw 


773x381x27x17 


754x268x17x13 


La carga correspondiente a los primeros 5 pisos es de 3.269 T/m y para el 
ultimo piso se asigna 1.852 T/m como se muestra en la figura 5. 
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Figura 5 Distribucion de cargas asignadas en las vigas del portico en analisis. 


3. DATOS PARA CEINCI-LAB 

3.1 Geometria del portico principal 

Con dos vectores se define la geometria de un portico regular en 
elevacion, en el sistema de computacion CEINCI-LAB. Estos vectores son: svque 
contiene la separacion de los vanos y sp la altura de cada uno de los pisos. Aguiar 
etal. (2019). 

A continuacion, en la figura 6 se muestra un fragmento del codigo del 
programa, donde la primera funcion que se utiliza es global seguido de una lista 
de variables, que son comunes para CEINCI-LAB y para OpenSees, se definen 
una sola vez, es similar a la funcion Common que se tiene en Fortran. 

Luego se muestra los datos de geometria del portico (sv , sp), pesos y 
cargas distribuidas en las vigas de cada piso (gp). Adicionalmente se emplean 
funciones conocidas de CEINCI-LAB para determinar las coordenadas de los 
nudos (X, Y ), la identificacion del nudo inicial y final de cada elemento (Nl, NJ) y 
los parametros como longitud, seno y coseno de cada barra que corresponden a 
columnas y vigas en el portico de analisis (L, seno, coseno). Tambien es 
importante conocer los grados de libertad que se deben asignar a cada nudo 
denominados CG y el vector de colocacion l/Cque es la asignacion de los grados 
de libertad en los nudos iniciales y finales de cada elemento. El marco teorico se 
presenta en Aguiar et al. (2019). 




















Acoplamiento de OpenSees con CEINCI-LAB ... 


375 


Con esto se tienen los datos con los cuales se calcula la matriz de rigidez 
del Portico, por ensamblaje directo y el vector de cargas generalizadas. 


Codigo del programa principal, parte 1: 


global DOF nod sv sp nv np nr X Y Nl NJ n masa_nudo Xcol Ypisos L coseno Es 
Abeam Ibeam Mybeam Ks_beamxL qp Acol lcol Mycol Ks_colxL num elem 
Tabla_Elem_Nud DP DN th pP th pN th pcP th pcN Dis_Fuerzas_Laterales T W fi 

%% Geometria del portico 

sv =[7.315;7.315;7.315]; %lngresar un vector con la longitud de vanos en 
(metros) 

sp =[5.486;3.81; 3.81; 3.81; 3.81; 3.81]; % Vector con la altura de pisos 
(metros) 

g=9.81; % Gravedad m/seg2 

qp=[-3.269;-3.269;-3.269;-3.269;-3.269;-1.8517]; % carga distribuida en vigas 
en cada piso T/m 

[nv,np,nudt,nudcol,nudvg,nod,nr]=geometria_nudo_viga(sv,sp); 

[X,Y]=glinea_portico2(nv,np,sv,sp,nod,nr); 

[NI,NJ]=gn_portico2(nr, nv, nudt, nudcol, nudvg); 
[CG,ngl]=cg_sismo2(nod,nr,Y); %Calcula los grados de libertad 
[VC]=vc(NI,NJ,CG); 

[L,seno,coseno]=longitud(X,Y,NI,NJ); %long. de @elem, seno y coseno 
Xcol=X(1,1 :nr); 

Ypisos=nonzeros(unique(Y)); 


Figura 6 Codigo del programa principal, parte 1. 

3.2 Definicion de masa 

En el analisis estructural la matriz de masa es importante, porque la 
respuesta del analisis depende de este parametro. Entre las opciones para 
modelar estas matrices se tiene, las matrices de masas concentradas (Lumped 
Mass Matrix) y matrices de masas consistentes (Consistent Mass Matrix). 

En una matriz de masas concentradas, la masa de una parte de la 
estructura analizada es asignada al grado de libertad asociado, que esta definido 
por un termino de la diagonal de la matriz (Garcfa, 2013). Para definir esta matriz 
existen varios metodos, entre ellos los propuestos en trabajos de: Zaradnik, 
(2009); Zaradnik et al. (2006); Hughes (1987); tambien es posible definirla 
directamente en algunos casos. 

Las matrices de masas consistentes, resultan de analisis variacional y 
tienen valores diferentes de cero, en la diagonal principal y fuera de ella (Paz, 
1992). Para la definicion de esta matriz se debe recurrir a la formulacion de 
elementos finitos como se presenta en la ecuacion 1 (Felippa, 2004), 
considerando la energfa cinetica del sistema e interpolando las velocidades, de 
forma similar a la que se interpolan los desplazamientos para la formulacion de la 
matriz de rigidez. 
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niij = j Ni T -p-NjdV (1) 

v 

Donde Ni y Nj corresponden a las funciones de forma elementales que se 
hayan adoptado; p es la densidad del material y V el volumen del elemento. 

Desde un punto de vista energetico, esta matriz de masa es consistente 
con la matriz de rigidez porque ambos se basan en las mismas funciones de 
desplazamiento (Carr, 2009). 

Acorde a Zaradnik et al. (2009), la definicion de la matriz de masa 
empleada en un modelo tiene un gran impacto en el esfuerzo computacional al 
momento de resolver problemas numericos. Ademas de Zaradnik, otros autores 
apoyan la premisa de que el modelo con masas consistentes podrfa llevar a 
resultados con mayor precision, pero el esfuerzo computacional serfa muy grande 
respecto al modelo donde se utilice masas concentradas. El modelo de masas 
concentradas ademas presenta la ventaja de que permite desacoplar las 
ecuaciones de movimiento de forma sencilla. Chopra (2012), indica que es mas 
empleado el modelo de masas concentradas. Entre los parametros de rigidez, 
masa, amortiguamiento para la simulacion estructural, Carr (2009) manifiesta que 
la representacion de rigidez es la mas significativa en el analisis, el modelo de 
masa no es tan critico siempre que se represente la masa total de la estructura. 

Para los modelos planteados en esta investigacion se emplea la matriz de 
masa concentrada (Lumped Mass Matrix), es decir que se asigna la masa de los 
elementos y cargas distribuidas a los nudos, de la union viga-columna, para 
facilitar el analisis con OpenSees, a continuacion, se muestra en la figura 7 la 
distribucion de las masas nodales. 


De forma general se calcula la carga colaborante a cada nudo (union viga- 
columna), multiplicando la carga uniforme distribuida por el ancho colaborante, 
para este caso se define, de forma simplificada, como la mitad de la luz de cada 
vano que se une al nudo considerado; luego se halla la masa, dividiendo esta 
carga para el valor de la aceleracion de la gravedad. Sumando los pesos de cada 
piso se halla el peso total del edificio. 
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Figura 7 Distribucion de masas en los nudos, para las funciones de OpenSees. 
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3.3 Patron de fuerzas laterales para Pushover monotonico 

El procedimiento para el analisis estatico no lineal, conocido como 
Pushover, se describe en el documento ASCE/SEI 41-13, seccion 7.4.3. De forma 
general, la aproximacion estatica consiste en aplicar fuerzas laterales a un modelo, 
considerando los efectos de no linealidad, estas fuerzas van aumentando hasta 
que se obtenga un desplazamiento objetivo o se alcance el colapso (Arango et al., 
2013). La curva de capacidad de una estructura se obtiene graficando los valores 
maximos de cortante basal y desplazamiento en el ultimo nivel, como se muestra 
en la figura 8. 




Figura 8 Secuencia de analisis Pushover. 


Segun Carr (2009), la capacidad de la estructura depende del patron de 
carga utilizado, es por ello que este es uno de los aspectos centrales en el analisis 
Pushover. La distribucion de fuerza sfsmica tfpica, que se encuentran en los 
codigos de disefio o que recomiendan autores en artfculos de investigacion, 
dependera de la masa en cada nivel. Esta masa se distribuira en la altura de la 
estructura en base a la siguiente ecuacion. 


w x h x 


ZU^tht 


( 2 ) 


Donde Cvx es un factor de distribucion vertical de fuerzas laterales; w t es 
la parte del peso total en la posicion i, es decir en el nivel i; w x es la porcion del 
peso total en la posicion x, de igual modo esta posicion representa el nivel x; es 
la altura del nivel i desde la base; h x es la altura del nivel x medido desde la base 
y el factor k depende del metodo de distribucion de fuerzas laterales. 


Para el modelo de analisis en este artfculo se debe definir como se 
realizara la distribucion de fuerzas laterales en cada piso de la edificacion, esto se 
lo indica con la variable Dis_Fuerzas_Laterales. El valor que esta variable puede 
tomar se encuentra entre 1 y 4, donde 1 corresponde a una distribucion en base al 
primer modo de vibracion de la estructura (en la ecuacion 2, k = 1 y en lugar de 
emplear el peso de la estructura se emplea el patron de deformada del primer 
modo de vibrar); 2 es una distribucion de fuerzas parabolica, donde k varia en 
funcion del periodo fundamental de la estructura, es decir que k = 1 para T < 0.5 
(s) o k = 2 para T > 2.5 (s), para periodos en el rango de 0.5 a 2.5 (s) se empleara 
una variacion lineal para estimar k en funcion del periodo fundamental; para un 
valor de 3, se tiene una distribucion triangular, en este caso k=1; finalmente si la 
variable Dis_Fuerzas_Laterales es 4 la distribucion sera uniforme, es decir que 
las fuerzas en cada piso seran iguales. 
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Para el analisis de la estructura, mostrada en la figura 4, se presenta la 
distribucion de fuerzas laterales en funcion de la altura de la edificacion, como se 
ilustra en la figura 9, mas adelante se presenta el calculo de su periodo 
fundamental que es igual a 0.7 (s). 



Figura 9 Distribucion de fuerzas laterales en funcion de la altura de la estructura. 

En la figura 10 se presenta el codigo con los datos de masas 
determinados en funcion de la carga distribuida en cada piso. Adicionalmente se 
escoge Dis_Fuerzas_Laterales=1, es decir, la distribucion de fuerza laterales del 
patron de cargas sera en base a la deformada del primer modo de vibracion. 


Codigo del programa principal, parte 2: 

%% Masas 
W=zeros(np,1); 
fort=1:np 
sum=0; 
for r=1 :nv 

sum=-qp(t,1)*sv(r,1); 

W(t,1)=W(t,1)+sum; 

end 

end 

sum=0; WBuilding=0; 
for i=1 :np 

WBuilding=WBuilding+W(i,1); %Peso total en T 
for ij=1 :nv+1 
if ij==1 

NodalMass =(-qp(i,1 )/g)*sv(ij,1 )/2; 
elseif ij==nv+1 

NodalMass =(-qp(i,1 )/g)*sv(nv,1 )/2; 
else 

NodalMass =(-qp(i,1 )/g)*sv(ij-1,1 )/2+(-qp(i,1 )/g)*sv(ij,1 )/2; 
end 

masa_nudo(i,ij)=NodalMass; 

end 

MASA(i,i)=W(i,1)/g; 

end 

Dis_Fuerzas_Laterales = 1; % 1 En funcion del primer modo de vibrar 


Figura 10 Codigo del programa principal, parte 2. 
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Con esta informacion es posible generar graficos del portico en analisis 
con la funcion dibujo, tambien se puede representar los elementos y nudos con su 
identificacion con dibujoNudElem como se muestra en la figura 11 b (notar que las 
vigas consideran un nudo en el medio, es decir que en cada vano del portico se 
modelan 2 vigas). La funcion dibujogdl genera un grafico con los grados de 
libertad que se han asignado a cada nudo como se presenta en la figura 11 c. 


Codigo del programa principal, parte 3: 


% Dibujo de la Estructura 2D 
dibujo(X,Y,NI,NJ) 
dibujoNudElem(X,Y,NI,NJ,CG) 
dibujogdl(X,Y,NI,NJ,CG) 
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Figura 11 a) Codigo para dibujar las estructuras; b) Identificacion de nudos y 
elementos; c) Grados de libertad para el analisis en CEINCI-LAB 


4. MODELO CON ZONA DE PLASTICIDAD CONCENTRADA 
DE OPENSEES 

Para el analisis Pushover no lineal en OpenSees de este artfculo se 
trabaja con un modelo de plasticidad concentrada, por ello es necesario generar 
un modelo con la ubicacion de las zonas de plasticidad concentrada asignando 
las propiedades de no linealidad de cada seccion. Ibarra etal. (2005) 

Para explicar este modelo se muestra un portico de 2 vanos y 1 piso, 
donde se tiene 3 columnas y 4 vigas, considerando que en cada vano existen 2 
vigas de igual longitud. La generacion del nudo en el medio del vano es util para 
incluir elementos diagonales de acero o disipadores de energfa como TADAS o 
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ADAS. Luego cada columna y cada viga mostrada en la figura 12a es dividida en 5 
partes como se muestra en la siguiente figura 12b donde a los extremos se tiene: 
elementos rigidos luego una rotula plastica de longitud nula seguido de un 
elemento elastico. Lignos y Krawinkler (2007, 2009, 2011), Lignos et al. (2011). 
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Figura 12 a) Modelo generado en CEINCI-LAB b) Modelo para el analisis 
Pushover no lineal en OpenSees. 

En el modelo mostrado anteriormente se tiene 4 materiales de los 
elementos vigas-columnas y son: elemento rigido; zona de dano o de plasticidad 
concentrada para las columnas de acero; zona de dano para viga; elemento 
elastico. En este analisis Pushover desarrollado en OpenSees se considera 
conveniente asignar en las zonas de plasticidad concentrada modelos bilineales, 
para columnas de acero con el material denominado SteelOI y para vigas se 
emplea el modelo denominado Ibarra-Medina-Krawinkler, existe un tramo del 
elemento que presentara un comportamiento elastico con el material Elastic beam 
column. Estos materiales se estudiaran mas adelante y se revisara los factores de 
modificacion de rigidez para el analisis de la estructura. 

En la figura 13 se puede visualizar la ubicacion de las zonas de plasticidad 
concentrada. Las zonas de plasticidad concentrada de las vigas se colocaron en 
las caras de la columna. Es decir, a una distancia de la mitad del ancho de las 
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columnas desde la Ifnea central de la columna. Ademas, se colocaron las zonas 
de plasticidad concentrada en columnas a una distancia de la mitad del ancho de 
la viga desde la linea central de la viga, como se lo realiza en el reporte de Shoma 
y Constantinou (2017). 

dcol dccH 



(a) (b) 

Figura 13 a) Zonas de plasticidad concentrada en vigas y columnas e 
identificacion de la zona rigida del elemento; b) Elementos elasticos. 

4.1 Elemento elastico 

Las columnas y vigas se modelan como elementos elasticos (Elastic 
beam column), en la zona central, donde los esfuerzos son directamente 
proporcionales a las deformaciones, como se muestra en la figura 14. 



Figura 14 Modelo de elementos elastico (comportamiento lineal). 

4.2 Modelo bilineal para zonas de plasticidad concentrada en 
columnas 

El modelo para el comportamiento no lineal de columnas es sencillo de 
entender, es un modelo bilineal con una rigidez elastica hasta un esfuerzo o fuerza 
denominado de fluencia y una zona de post fluencia posterior a ese nivel de 
solicitacion; este modelo generalmente se lo conoce como elastoplastico y en 
OpenSees se puede modelar con el comando steelOI. 
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Figura 15 Modelo Elastoplastico para las articulaciones que se pueden formar en 
la columna. 


4.3 Modelo de Ibarra-Medina-Krawinkler para zonas de plasticidad 
concentrada en vigas 

Este modelo es mas complejo que los indicados anteriormente, ya que 
tiene relacion con la degradacion de resistencia del material, la curva constitutiva 
del modelo se indica en la figura 16. En este artfculo se realiza un analisis 
Pushover monotonico, es decir que no existe un efecto reversible de la carga y 
por ello el modelo Ibarra-Medina-Krawinkler desprecia el efecto ciclico en el 
diagrama momento - rotacion, como la Ifnea en azul de la figura 16. 

Momento 



Figura 16 Modelo Momento-Rotacion de Ibarra-Medina-Krawinkler. 

En CEINCI-LAB, se han presentado modelos constitutivos de acero para 
zonas de plasticidad concentrada en base a las recomendaciones del ASCE/SEI 
41. En Aguiar et al. (2015), Mora y Aguiar (2015, 1), se describe el procedimiento 
de calculo para el diagrama de momento-rotacion de la seccion W 30x99, en este 
caso se presentaran las ecuaciones necesarias para desarrollar este diagrama en 
base al modelo de Ibarra-Medina-Krawinkler, que se abrevia como IMK. 
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Las expresiones necesarias para determinar los parametros del modelo 
momento-rotacion IMK dependen de la conexion que existe en las vigas con las 
columnas; de las conexiones precalificadas en el AISC 358 se tiene: 






(g) 

Figura 17 (a) Conexion de Patin Soldado Sin Refuerzo y Alma Soldada (WUF-W); 

(b) Conexion de Seccion de Viga Reducida (RBS); (c) Conexion de 
Placa Extrema Empernada Sin Rigidizar (BUEP) y (d) rigidizada 
(BSEP); (e) Conexion de Placa Empernada (BFP); (f) Conexion Kaiser 
con Soporte Atornillado y (g) Soldado (KBB). 

En las siguientes ecuaciones se requiere ingresar el valor de Fy en 
unidades de Megapascal y para parametros de longitud debe utilizarse milimetros, 
en estos casos los coeficientes Cunitl y Cunit2 son 1. La relacion de Mc / My es, en 
promedio, 1.10. 
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■ Conexion viga-columna diferente a la RBS: 

La rotacion que existe en el rango de fluencia se determina en base a la 
ecuacion 3 y la rotacion en el rango de post fluencia con la ecuacion 4. 



Donde 0p, es la rotacion plastica pre-capping (rad); 0pc, es la rotacion 
plastica post-capping (rad); d, altura total de la viga (mm); h, es la altura del alma 
de la viga (mm); tw, es el espesor del alma de la viga (mm); bt, es el ancho del 
patin de la viga (mm); tf, es el espesor del patin de la viga (mm); L, es la longitud 
de la viga (mm); Fy, es el esfuerzo de fluencia (MPa); Cunit^l ; Cunit 2 =1. 

■ Conexion viga-columna RBS: 

De manera similar al caso de vigas con conexion diferente a la RBS, se 
tiene la ecuacion 5 y 6 para determinar la rotacion que existe en el rango de 
fluencia y de post fluencia respectivamente. 



Algunos terminos de las ecuaciones fueron descritos previamente, 
restando indicar: L b , que es la longitud de la viga sin arriostramiento (mm); r y es el 
radio de giro en (mm). Por otra parte, se debe cumplir: 

L„ 

~ =£ 70 (?) 

'y 

Para el modelo IMK en vigas se puede entender de manera conservadora 
que el valor de My es el producto del modulo plastico de la seccion por el esfuerzo 
de fluencia del material, sin embargo, cuando se emplean conexiones 
precalificadas este valor puede ser multiplicado por un factor que toma en cuenta 
el endurecimiento por deformacion denominado Cpr, siendo este valor calculado 
con la siguiente ecuacion. 


F + F 

Cpr = |- — - ; Cpr < 1.2 ( 8 ) 

Donde Fy, es el esfuerzo de fluencia del material (ksi); Fu, es el esfuerzo 
ultimo del material (ksi). 
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El valor de la rotacion donde termina el rango elastico se denomina 9 y , 
este valor puede ser calculado en base a la ecuacion del ASCE/SEI 41, como: 


M„ ■ L 

e y = p — 

V 6 ■ E ■ Iy 


(9) 


Donde, M p es el momento plastico (T.m); L es la longitud de la viga en (m); 
E, es el modulo de elasticidad del material; / y , es el momento de inercia de la 
seccion respecto al eje de giro (m 4 ). 


Para calcular el momento residual, M r , 
de la ecuacion 10. 


Mr 

My 


0.4 


se recomienda utilizar la relacion 

( 10 ) 


Se genera un programa para comparar los modelos de momento rotacion 
de IMK con el ASCE/SEI 41, este programa forma parte de CEINCI-LAB y se 
denomina Momento Rotacion Curvtura ASCE41 Acero BCG.m 


Un fragmento de los datos de este programa se presenta en la figura 18, 
en esta parte se requiere datos del material como, modulo de elasticidad, esfuerzo 
de fluencia del acero y esfuerzo ultimo, ademas, de la longitud de la viga y la 
estimacion de la longitud plastica de la seccion. 


Momento Rotacion Curvtura ASCE41 Acero BCG.m 


% Materiales 

Ey=200.0*1000000.0/9.81 ;% Modulo de elasticidad, T/m2 
fy=25310.505; % Esfuerzo de fluencia, T/m2 
fu=40778.036; % Esfuerzo ultimo, T/m2 
L = 7.315; % Longitud, m 
Lp=0.15; % Longitud plastica, m 

%.W30x99 Viga 

d = 754.38; % Peralte, mm 

tw = 13.208; % Espesor del alma, mm 

bs = 267.7; % Patfn superior, mm 

bi = 267.7; % Patfn inferior, mm 

tfs = 17.018; % Espesor patfn superior, mm 

tfi = 17.018; % Espesor patfn inferior, mm 


Figura 18 Datos para Momento - Rotacion modelos ASCE/SEI 41 y IMK. 

Para una seccion W30x99, los modelos momento rotacion se presenta 
en la figura 19. En este caso es evidente que los modelos del ASCE/SEI 41 
difieren del modelo IMK, siendo mas conservador el modelo del ASCE/SEI 41. 
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Diagrama Momento * Rotacibn Diagrama Momento - Rotacion 

tHUr---»-»------ -1 140 (------r-i- 1 - 1 - 



Figura 19 a) Momento-Rotacion con conexion RBS; b) Momento-Rotacion con 
conexion diferente que RBS. 

Para el portico en este modelo se asume que la conexion de la viga a la 
columna es RBS, los valores de las variables “DP, DN, th pP, th_PN, th_pcP, 
th pcN” se presentaran mas adelante para las zonas de plasticidad concentrada 
en las vigas con el modelo de Ibarra-Medina-Krawinkler. 

Donde DP, es tasa de deterioro ciclico en la direccion de carga positiva; 
DN, es tasa de deterioro ciclico en la direccion de carga negativa; DP y DN son 
parametros que se utilizan para crear un comportamiento histeretico asimetrico; 
th pP = 0p, es la rotacion en el rango de fluencia en el caso positivo (rad); th_pN = 
0p, es la rotacion en el rango de fluencia en el caso negativo (rad); th_pcP = 0pc, 
es la rotacion en el rango de post fluencia en el caso positivo (rad); th pcN = 0pc, 
es la rotacion en el rango de post fluencia en el caso negativo (rad). 


4.4 Modificacion de rigidez de los elementos en el analisis con 
OpenSees 

En este modelo de plasticidad concentrada cada elemento, viga o 
columna, es considerado como un elemento elastico, conectado en serie por 
zonas de plasticidad concentrada como se muestra en la figura 20. 

Ks Kbc Ks 


Figura 20 Distribucion de rigidez en un elemento. 

Acorde a la tesis de Barrera (2018), las propiedades estructurales de cada 
miembro son una combinacion de propiedades de tres subelementos; estos 
elementos son, las zonas de plasticidad concentrada en los extremos con 
propiedades no lineales y del elemento viga o columna con propiedades elasticas. 
Sin embargo, para mantener la rigidez equivalente del conjunto de forma 
consistente con la rigidez del elemento verdadero es necesario modificar la rigidez 
de estos subelementos. 
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El modelo idealizado en la figura 20, muestra que la rigidez de cada 
subelemento se puede relacionar con la rigidez de un elemento de la siguiente 
manera: 


1 _ 1 1 
Kelem Ks Kbc 


( 11 ) 


Ademas, la rigidez rotacional de un elemento sometido a doble curvatura 
se puede determinar por: 


Kelem = 


6 ■ E ■ l y 
Lelem 


( 12 ) 


Con el proposito de resolver la ecuacion de rigidez en serie de los 
subelementos que conforman a cada miembro y para evitar problemas de 
inestabilidad inelastica numerica, los autores Ibarra y Krawinkler proponen usar 
modelos con zona de plasticidad concentrada con una rigidez “n” veces mayor que 
la rigidez rotacional del elemento elastico (Ibarra & Krawinkler, 2005), 
generalmente se emplea valores de n cercanos a 10. 

De esta forma, la rigidez elastica de las zonas de plasticidad concentrada 
(Ks) y del elemento elastico (Kbc) deben ser modificadas por este valor n, 
considerando las siguientes ecuaciones: 

n +1 .... 

K bc - - Kelem (13) 

n 

K s = (n + i)- Kelem (14) 

5. DATOS PARA OPENSEES GENERADOS POR CEINCI-LAB 

Se debe ingresar datos del material y secciones para el Pushover, se 
asigna las propiedades como se indica en el siguiente fragmento del programa: 


Codigo del programa principal, parte 4: 

%% Propiedades del material 
Fy=25310.505; % Esfuerzo de fluencia, T/m2 
Es=200.0*1000000.0/9.81;% Modulo de elasticidad, T/m2 
Fu=40778.036; % Esfuerzo ultimo, T/m2 

Cpr=(Fy+Fu)/(2*Fy); % Factor que toma en cuenta la resistencia maxima de la 
conexion. 
if Cpr < 1.2 
Cpr=Cpr; 
else 

Cpr =1.2; 

End 

%% Secciones de los elementos, Columna % W30x173 
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Zx = 0.0098; 

Acol = 0.0326; 
lcol = 0.0034; 

Mycol =Fy*Zx*Cpr; 
d_col= 773.2/1000.0; 

% Viga % W30x99 

Zx = 0.0051; 

Abeam = 0.0186; 

Ibeam = 0.0017; 

Mybeam = Fy*Zx*Cpr; 
d_beam = 754.4/1000.0; 

% Valores para zonas de dano del Modelo IMK 
DP= 0.886; % rate of cyclic deterioration for pos loading 

DN= 0.886; % rate of cyclic deterioration for neg loading 

th_pP= 0.0204; % plastic rot capacity for pos loading 
th_pN= 0.0204; % plastic rot capacity for neg loading 
th_pcP= 0.1001; % post-capping rot capacity for pos loading 
th_pcN= 0.1001; % post-capping rot capacity for neg loading 


Figura 21 Codigo del programa principal, parte 4. 

Por facilidad en este primer programa, las vigas y columnas son de 
seccion unica, la seccion de columnas es W30x173 y las vigas W30x99, para 
cumplir con este proposito se emplea la funcion denominada 
[Seccion]=Unica_Vg_COL(nudt, coseno, Acol, lcol, Abeam, Ibeam); los datos 
que requiere son numero total de vigas mas columnas nudt, el area de la seccion 
transversal e inercia de la seccion de columnas Acol e lcol, de igual modo 
respecto a la seccion de las vigas Abeam, Ibeam y con ello se genera una matriz 
denominada Seccion. La matriz Seccion ya se manejaba anteriormente en 
programas de CEINCI-LAB. Adicionalmente se determina la matriz de rigidez total 
de la estructura, el periodo de vibracion y propiedades dinamicas de la estructura 
como modos de vibracion y factores de participacion modal. 


Codigo del programa principal, parte 5: 


%% Propiedades dinamicas 

[Seccion]=Unica_Vg_COL(nudt,coseno,Acol,lcol,Abeam,lbeam); 

[ELEM]=gelem_portico(Seccion); 

[K]=krigidez_acero(ngl,ELEM,L,seno,coseno,VC,Es); 

KTT=K; %<.MATRIZ "K" DE TODA LA ESTRUCTURA (T/m) 

na=np; 

kaa=KTT(1 :na,1 :na);kab=KTT(1 :na,na+1 :ngl); 
kba=kab'; kbb=KTT(na+1 :ngl,na+1 :ngl); 

KL=kaa-kab*inv(kbb)*kba; 

[T,fi,OM]=orden_eig(KL,MASA); 

% Factores de participacion modal 

b=ones(np,1); 

for i=1 :np 
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maxfi=max(abs(fi(:,i))); 

fi(:,i)=fi(:,i)/maxfi; 

end 

for i=1 :np 

FP(:,i)=abs(((fi(:,i))'*MASA*b)/((fi(:,i))'*MASA*(fi(:,i)))); 

end 


Figura 22 Codigo del programa principal, parte 5. 

Regresando al portico de analisis de 6 pisos mostrado en la figura 4, es 
necesario generar un modelo con la division de los elementos en las partes 
correspondientes como se indica en la figura 12. 


Para el modelo de OpenSees se requiere nuevos nudos, coordenadas X, 
Y e identificacion de los nudos iniciales y finales que componen a cada elemento, 
a mas de ello es necesario identificar las secciones que se deben asignar a cada 
elemento. Se genera la funcion [Tabla_Elem_Nud, nod, X, 
Y]=ElementosNudos(nudt, nudcol, nod, Nl, NJ, X, Y, d_col, d_beam); en este 
caso la matriz Tabla_Elem_Nud contiene la informacion relevante de nudos y 
elementos. Generando un modelo como se indica en la figura 23. 
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(a) 


(b) 


Figura 23 a) Numeracion de elementos columnas; b) Numeracion de elementos 
vigas. 


En este nuevo modelo se puede visualizar 300 elementos, es por ello que 
la matriz denominada Tabla_Elem_Nud tiene 300 filas donde cada una de ellas 
corresponde a un elemento. Esta matriz tiene 4 columnas donde la primera es una 
etiqueta, es decir el primer numero corresponde al elemento original (figura 11 a) 
seguido de 2 ceros y el ultimo numero, del 1 al 5, representa el segmento del 
elemento que fue dividido (figura 12b); la columna 2 y 3 corresponde al nudo inicial 
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y final de cada nuevo elemento, la columna 4 es una etiqueta para identificar el 
tipo de propiedad que corresponde a cada elemento. La etiqueta de la columna 4 
puede ser 0 para elementos rigidos (misma propiedad que un elemento elastico, 
pero con area e inercia multiplicada por 10), luego 1 para elementos elasticos, 2 
para zonas de plasticidad concentrada con propiedades de un elemento bilineal 
steelOI y 3 para las zonas de plasticidad concentrada de vigas con propiedad del 
modelo Ibarra-Medina-Krawinkler. 
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(b) 


Figura 24 Fragmento de Tabla_Elem_Nud a) Primeras 2 columnas (elementos 1 y 
2 de figura 11 a) b) Primeras 2 vigas (elementos 25 y 26 de figura 11a). 

Un fragmento del programa se detalla en la figura 25, en el cual se puede 
notar una nueva funcion denominada Pushover_OpenSees(), que se presenta de 
color rojo para distinguirla, en este caso no existen datos de ingreso y salida 
debido a que esta rutina trabaja con variables globales definidas al inicio del 
programa. Sin embargo, al ejecutarse esta linea se genera un archivo con 
extension TCL cuyo nombre es Estructura_Pushover.tcl. 


Codigo del programa principal, parte 6: 


%% Elementos y nudos 

[Tabla_Elem_Nud,nod,X,Y]=ElementosNudos(nudt,nudcol,nod,NI,NJ,X,Y,d_col,d. 

beam); 

NI=Tabla_Elem_Nud(:,2)'; 

NJ=Tabla_Elem_Nud(:,3)'; 

[L,seno,coseno]=longitud(X,Y,NI,NJ); %long. de @elem, seno y coseno 
num_elem=nudt*5; 

HBuilding=max(Y); 

%% Nudos con grados de libertad iguales 
[DOF]=MismosGDL(nod,np,Xcol,Ypisos,X,Y); 

Pushover_OpenSees() 


Figura 25 Codigo del programa principal, parte 6. 
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En la figura 26 se puede visualizar el funcionamiento de la funcion 
Pushover_OpenSees(), los datos de ingreso corresponden al codigo generado 
hasta este punto en CEINCI-LAB y como salida se tiene un archivo con el modelo 
en lenguaje TCL. Mas adelante se presentara un detalle de los modelos de 
transformacion geometrica de los elementos estructurales empleados. 



iwo TCL con » wnlHo 

O Ort nx i p mm i on rn n 



Figura 26 Diagrama de proceso de la funcion Pushover_OpenSees() 

Como una forma de verificar graficamente el modelo generado hasta este 
punto se puede elaborar los dibujos del portico en analisis, en este caso la funcion 
dibujoNudElemMODPUSH(X,Y,NI,NJ) presenta un grafico como se muestra en la 
figura 27b, en el cual se incluyen los nudos y elementos columnas y vigas que 
componen al portico, con una distribucion similar a la presentada en la figura 23. 


Codigo del programa principal, parte 7: 


%% Dibujos 

dibujoNudElemMODPUSH(X,Y,NI,NJ) 


(a) 


Ciqurcmj Eslructural 



(b) 


Figura 27 a) Codigo del programa principal, parte 7. b) Modelo de vigas y 
columnas generado en CEINCI-LAB. 
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El archivo Estructura_Pushover.tcl que se genera en la carpeta que se 
esta almacenando el programa en MATLAB, es ejecutado posteriormente por 
OpenSees y genera los calculos necesarios para el Pushover, la linea que indica 
esta accion se muestra a continuacion, donde el comando empieza con un signo 
de admiracion (!) luego se escribe OpenSees.exe seguido por el nombre del 
archivo tcl que describe el modelo, en este caso Pushover_Estructura.tcl 


Codigo del programa principal, parte 8: 


%% Analisis Pushover 

! OpenSees.exe Pushover_Estructura.tcl 


Figura 28 Codigo del programa principal, parte 8. 

Finalmente, al ejecutarse OpenSees desde Matlab se generan los calculos 
del Pushover del modelo propuesto y aparece una ventana emergente con la 
simulacion del analisis como se muestra en la figura 29. 



Figura 29 Ventana con la simulacion del analisis Pushover en OpenSees. 

Para que el analisis se desarrolle de forma correcta, en la carpeta junto al 
programa Pushover_Estructura.tcl se debe encontrar los archivos con los 
modelos de materiales y tipos de elementos, a mas de los programas con 
extension (.tcl) necesarios para el analisis de OpenSees.exe, es decir los archivos 
que se muestran en la figura 30 (estos programas seran almacenados en los 
perfiles de ResearchGate de los autores cuando se publique el artfculo). 


# •'usrwrti.euruoura Programa generado desde CEINCI-LAB con el modelo 


# DispljyPUm* 

# Mmom 

# prcntnode 

# roltrjningCoi 

# RotSpnng2DMOdlKMoO*l 

# RotSpnngSlrdOI 

# RotSpnng$t««l02 


Programas necesarios para que se ejecute el analisis. 


OprnSw 


Programa OpenSees.exe 


Figura 30 Programas necesarios para que se ejecute el analisis en OpenSees. 
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Los resultados del analisis Pushover se almacenan en una carpeta que es 
generada de forma automatica por OpenSees, esta carpeta tiene por nombre 
Concentrated-Pushover-Output y almacena en este caso 5 archivos de texto con 
extension .out, estos archivos son Desplazamiento en el tope Disp_Roof, y para 
cada nudo de la base se genera un archivo denominado Vbase# con informacion 
del cortante experimentado como reaccion ante las cargas laterales aplicadas. 


La funcion Graficar_Pushover2.m es la que genera la grafica con 
desplazamiento en el tope y cortante basal, es decir la curva del analisis Pushover, 
como se muestra en la figura 31 b. La metodologfa para determinar el 
desplazamiento efectivo de fluencia S y y desplazamiento ultimo S u es la que se 

presenta en el FEMA P695, donde el desplazamiento ultimo es el valor de 
desplazamiento para el cual se ha tenido una perdida del 20% de resistencia, es 
decir el valor para el cual el cortante baja a 0.8 ■ Vmax. El desplazamiento efectivo 
de fluencia se puede calcular en base a la ecuacion 15. 


y eff 


= n ■ 


^max / 9 \ 
W 'U■ n 2 ' 


■T^ 


(15) 


Donde n es el factor de participacion modal del primer modo de vibracion 
bajo condiciones elasticas; T x es el periodo fundamental de vibracion de la 
estructura; g es la gravedad; V mkx es el cortante maximo del analisis Pushover; W 
es la carga reactiva total del sismo. 


Codigo del programa principal, parte 9: 

%% Grafica Pushover 

Tfundamental=textread('Resultados_Modal\Periodo_Fundamental.out'); 

[Vbase,DT]=Graficar_Pushover2(nr,WBuilding,HBuilding,FP,Tfundamental); 


(a) 


Curva Pushover 



(b) 

Figura 31 a) Codigo del programa principal, parte 9. Curva Pushover generada 
por CEINCI-LAB mediante el acoplamiento con OPENSEES 
(Distribucion de cargas en funcion del 1 er modo de vibracion). 
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6. REFORZAMIENTO DE LA ESTRUCTURA MEDIANTE 
DIAGONALES DE ACERO 

Muchas veces los porticos de acero, como el presentado en la figura 4, 
son reforzados con diagonales de acero con el proposito de aumentar la rigidez y 
resistencia ante solicitaciones sfsmicas. 

Las configuraciones de las diagonales de acero habituales en Ecuador 
para el reforzamiento de estructuras son las denominadas en forma de “V 
invertida” o “X de multiples pisos”. A continuacion, se presentara la metodologfa 
para realizar el analisis estatico no lineal (Pushover) de las estructuras con 
diagonales de acero presentadas en la figura 32. En el vano central se colocara 
diagonales de acero circulares HSS7x0.315, cuyas dimensiones son 18cm (aprox. 
7 pulgadas) de diametro exterior y 8mm (0.315 pulgadas) de espesor. 
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(b) 

Figura 32 a) Estructura con diagonales en V invertida; b) Estructura con 
diagonales en X de multiples pisos. 


7. DATOS PARA CEINCI-LAB INCLUYENDO DIAGONALES 

Al inicio del programa se debe declarar variables globales, en este caso 
existen nuevas variables, adicionales a las previamente ingresadas, entre ellas: 


Codigo del programa principal, parte 1b: 

global nod_Diag X_Diag Y_Diag NI_Diag NJ_Diag Abrace Ibrace num_Diag 
Tabla_Elem_Nud_Diag 

global Fy_brace Es_brace Dis_Fuerzas_Laterales T W fi 


Figura 33 Codigo del programa principal, parte 1b. 
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Para definir la geometrfa del portico conformada por columnas y vigas se 
presento una metodologfa que requerfa de los vectores sv y sp, para la 
separacion de vanos y la altura de cada uno de los pisos Aguiar et al. (2019). A 
continuacion, se debe incluir datos de los marcos en los cuales se va a colocar las 
diagonales en forma de V invertida y en forma de V para el modelo con X de 
multiples pisos. Por facilidad se debe numerar los marcos como se muestra en la 
figura 3, empezando en 1 desde la esquina inferior izquierda hacia la derecha, 
hasta el ultimo marco del piso 1, luego en los siguientes niveles se continua con la 
numeracion en el mismo sentido. Para la elaboracion de la figura 3 se puede 
emplear la funcion de CEINCI-LAB denominada dibujovanos(sv, sp) presentada 
en Herrera et al. (2019). 

Para la generacion del modelo con V invertidas de la figura 32a, se debe 
incluir el vector mar, este vector debe contener el numero del marco en cual se 
colocaran dichas diagonales. Posteriormente se debe definir el nudo inicial y final 
de los elementos diagonales con la rutina geometri'a_nudo_diagonales y 
combinar los vectores con la geometrfa de los elementos totales, las funciones 
diferentes al programa del modelo sin refuerzo se han marcado de color rojo para 
que sea facil identificarlas y aplicar en los programas de los usuarios. 


Codigo del programa principal, parte 2b: 

%% Geometrfa del portico 

sv =[7.315;7.315;7.315]; %lngresar un vectorcon la longitud de vanos en (metros) 
sp =[5.486;3.81; 3.81; 3.81; 3.81; 3.81]; %lngresar un vector con la altura de pisos 
en (metros) 
dibujovanos(sv,sp); 

% Diagonales de acero, ubicacion y propiedades de la seccion 
mar = [2;5;8;11;14;17]; 

%Para diagonales se calcula el area e inercia 
diae=0.18; %Diametro exterior 
diai=0.18-8*2/1000; %Diametro interior 
g=9.81; % Gravedad m/seg2 

qp=[-3.269;-3.269;-3.269;-3.269;-3.269;-1.8517]; % carga distribuida en vigas en 
cada piso T/m 

[nv,np,nudt,nudcol,nudvg,nod,nr]=geometria_nudo_viga(sv,sp); 

[X,Y]=glinea_portico2(nv,np,sv,sp,nod,nr); 

[NI,NJ]=gn_portico2(nr, nv, nudt, nudcol, nudvg); 

[CG,ngl]=cg_sismo2(nod,nr,Y); %Calcula los grados de libertad 
[VC]=vc(NI,NJ,CG); 

[L,seno,coseno]=longitud(X,Y,NI,NJ); %long. de @elem, seno y coseno 

num_COL_VG=length(L); 

Xcol=X(1,1 :nr); 

Ypisos=nonzeros(unique(Y)); 

%. 

[GEN]=geometria_nudo_diagonales(nv, np, nudt, mar); 
[NI2,NJ2]=gn_portico(GEN); 

[NI_Diag,NJ_Diag]=gn_portico3(NI,NJ, NI2,NJ2); 


Figura 34 Codigo del programa principal, parte 2b. 
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Para la generacion del modelo con X de multiples pisos de la figura 32b, 
se debe incluir los vectores mar2 y marv. Donde mar2, describe los marcos que 
tienen diagonales en forma de V invertida y marv los que tienen diagonales en 
forma de V para que al unir estos elementos se formen las X de varios pisos. En 
este caso la funcion a emplear es geometri'a2_nudo_diagonales, como se 
muestra a continuacion en un fragmento del programa. 


Codigo del programa principal, parte 2b: 


%% Geometrfa del portico 

sv =[7.315;7.315;7.315]; %lngresar un vector con la longitud de vanos en (metros) 

sp =[5.486;3.81; 3.81; 3.81; 3.81; 3.81]; %lngresar un vector con la altura de pisos 

en (metros) 

dibujovanos(sv,sp); 

g=9.81; % Gravedad m/seg2 

qp=[-3.269;-3.269;-3.269;-3.269;-3.269;-1.8517]; % carga distribuida en vigas en 
cada piso T/m 

%. 

mar2 = [2;8;14]; 
marv = [5;11;17]; 

%Para diagonales se calcula el area e inercia 
diae=0.18; %Diametro exterior 
diai=0.18-8*2/1000; %Diametro interior 

%. 

[nv,np,nudt,nudcol,nudvg,nod,nr]=geometria_nudo_viga(sv,sp); 

[X,Y]=glinea_portico2(nv,np,sv,sp,nod,nr); 

[NI,NJ]=gn_portico2(nr, nv, nudt, nudcol, nudvg); 

[CG,ngl]=cg_sismo2(nod,nr,Y); %Calcula los grados de libertad 
[VC]=vc(NI,NJ,CG); 

[L,seno,coseno]=longitud(X,Y,NI,NJ); %long. de @elem, seno y coseno 
num_COL_VG=length(L); 

Xcol=X(1,1 :nr); 

Ypisos=nonzeros(unique(Y)); 

%. 

[GEN]=geometria2_nudo_diagonales(nv, np, nudt,mar2,marv); 
[NI2,NJ2]=gn_portico(GEN); 

[NI_Diag,NJ_Diag]=gn_portico3(NI,NJ, NI2,NJ2); 


Figura 35 Codigo del programa principal, parte 2b. 

En estos modelos la aplicacion de masas y distribucion de cargas laterales 
se mantiene como en el modelo del portico sin reforzar, que se puede visualizar en 
la figura 4. Es posible numerar los elementos del modelo de CEINCI-LAB y los 
grados de libertad con las funciones dibujoNudElem y dibujogdl, para visualizar 
la geometria que se tiene hasta este punto con cada modelo (Figura 36), se debe 
indicar que para las funciones se requiere hacer referencia a NI_Diag y NJ_Diag 
que son los vectores que tienen los nudos inicial y final de cada elemento 
correspondientes a columnas, vigas y diagonales como se muestra en el siguiente 
fragmento de codigo de la figura 36a. 
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Codigo del programa principal, parte 3b: 


% Dibujos 

dibujoNudElem(X,Y,NI_Diag, NJ_Diag,CG) 
dibujogdl(X,Y,NI_Diag, NJ_Diag,CG) 
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Figura 36 Estructura con diagonales en V invertida: a) Codigo del programa 
principal, parte 3b. b) Numeracion de elementos; c) Grados de 
libertad. Estructura con diagonales en X de multiples pisos: d) 
Numeracion de elementos; e) Grados de libertad. 

Para las propiedades de los elementos diagonales (brace), se deben 
incluir el modulo de elasticidad, el esfuerzo de fluencia del material que los 
conforman y calcular el area e inercia de estos elementos como se indica a 
continuacion: 

































398 


Brian Cagua, Roberto Aguiar, Julia Pilatasig y David Mora 


Codigo del programa principal, parte 4b: 

[VC_Diag]=vc(NI_Diag,NJ_Diag,CG); 

[L_Diag,seno_Diag,coseno_Diag]=longitud(X,Y,NI_Diag,NJ_Diag); %long. de 

@elem, seno y coseno 

%% Propiedades del material 

Fy=25310.505; % Esfuerzo de fluencia, T/m2 

Es=200.0*1000000.0/9.81 ;% Modulo de elasticidad, T/m2 

Fy_brace=Fy; 

Es_brace=Es; 

Fu=40778.036; % Esfuerzo ultimo, T/m2 

Cpr=(Fy+Fu)/(2*Fy); % Factor que toma en cuenta la resistencia maxima de la 
conexion. 
if Cpr < 1.2 
Cpr=Cpr; 
else 

Cpr =1.2; 

End 

% Column section % W30x173 
Zx = 0.0098; 

Acol = 0.0326; 
lcol = 0.0034; 

Mycol =Fy*Zx*Cpr; 
d_col= 773.2/1000.0; 

% Brace section % Tcircular Diametro de 18 cm y espesor de 8 mm 

Abrace=pi*diae A 2/4-pi*diai A 2/4; % area 
lbrace=pi*diae A 4/64-pi*diai A 4/64; % inercia 
% beam section % W30x99 
Zx = 0.0051; 

Abeam = 0.0186; 

Ibeam = 0.0017; 

Mybeam = Fy*Zx*Cpr; 
d_beam = 754.4/1000.0; 

% beam springs at Floor 2 # % W30x99 
DP= 0.886; % rate of cyclic deterioration for pos loading 

DN= 0.886; % rate of cyclic deterioration for neg loading 

th_pP= 0.0204; % plastic rot capacity for pos loading 
th_pN= 0.0204; % plastic rot capacity for neg loading 
th_pcP= 0.1001; % post-capping rot capacity for pos loading 
th_pcN= 0.1001; % post-capping rot capacity for neg loading 


Figura 37 Codigo del programa principal, parte 4b. 

En CEINCI-LAB se puede calcular la matriz de rigidez lateral del portico 
con la contribucion de los elementos de vigas y columnas de acero (KH), luego se 
suma la contribucion de la rigidez de las diagonales de acero (KA). En el codigo 
presentado a continuacion, el color azul representa las Ifneas adicionales de 
programacion debido a las diagonales de acero y en rojo se muestra la suma de 
las matrices de rigidez. 
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Codigo del programa principal, parte 5b: 


%% Propiedades dinamicas 

[Seccion]=Unica_Vg_COL(nudt,coseno,Acol,lcol,Abeam,lbeam); 

Seccion=[Seccion; 

length(NI)+1 Abrace Ibracell 1]; 

[ELEM]=gelem jDortico(Seccion); 

% Matriz de Rigidez de columnas y vigas 
ELEMH=ELEM(1 :length(NI),:); 

[KH]=krigidez_acero(ngl,ELEMH,L,seno,coseno,VC,Es); 

%Contribucion de elementos Diagonales de Acero 
ELEMA=ELEM(length(NI)+1 :end.:): 

LA=L_Diag(length(NI)+1 :end); senoA=seno_Diag(length(NI)+1 :end); 
cosenoA=coseno_Diag(length(NI)+1 :end); VCA=VC_Diag(length(NI)+1 :end,:); 
num_Diag=length(LA); 

Est = Es/10; % Modulo de elasticidad, (kg/cm2) 

H=sp*100; 

hci=d_col*100;hcd=d_col*100;hdis=0;Lon=sv(2,1)*100;hv=d_beam*100; 

Kef=zeros(np,1); 

for ie=1 :np 

[Kdiag(ie),Keq(ie)]=rigidez_equivalente_TADAS(hci,hcd,hdis,Lon,hv,H(ie),Kef 

(ie),Abrace*10000,Est); 

end 

Kdiag=Kdiag/10; % Para tener en T/m 
ELEMT_Diag=[Kdiag(1);Kdiag(1); 

Kdiag(2);Kdiag(2); 

Kdiag(3);Kdiag(3); 

Kdiag(4);Kdiag(4); 

Kdiag(5);Kdiag(5); 

Kdiag(6);Kdiag(6)]; 

[KA]=krigidez_tadas(num_Diag,ngl,ELEMT_Diag,senoA,cosenoA,VCA); 

%Rigidez de Diagonales 

KTT=KH+KA; %<-MATRIZ "K" DE TODA LA ESTRUCTURA (T/m) 

na=np; 

kaa=KTT(1 :na,1 :na);kab=KTT(1 :na,na+1 :ngl); 
kba=kab'; kbb=KTT(na+1 :ngl,na+1 :ngl); 

KL=kaa-kab*inv(kbb)*kba; 

[T,fi,OM]=orden_eig(KL,MASA); 

% Factores de participacion modal 

b=ones(np,1); 

for i=1 :np 

maxfi=max(abs(fi(:,i))); 
fi(:,i)=fi(:,i)/m axfi; 
end 

for i=1 :np 

FP(:,i)=abs(((fi(:,i))'*MASA*b)/((fi(:,i)) , *MASA*(fi(:,i)))); 

End 

disp('Los periodos de vibracion en segundos con CEINCI-LAB son:') 

T % seg 


Figura 38 Codigo del programa principal, parte 5b. 
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El modelo de CEINCI-LAB para el analisis de un portico en dos 
dimensiones es un modelo simplificado de 1 grado de libertad (gdl) por piso, en las 
figuras 11 c, 36c y 36e se puede visualizar los grados de libertad considerados y se 
indica de color azul la componente principal de los grados de libertad horizontales. 
Para el calculo en CEINCI-LAB se debe realizar una condensacion estatica en los 
grados de libertad asociados a las componentes horizontales para generar el 
modelo simplificado. Este modelo para un portico de 6 pisos se presenta en la 
figura 39. 


m6 

• ► 6 

m5 


m4 

• -4 

m3 

•—► 3 

m2 

•—► 2 

ml 

• ► 1 


Figura 39 Modelo simplificado de 1 gdl considerado en CEINCI-LAB para el 
analisis. 


Para los 3 modelos del portico en analisis se ilustrara la matriz de rigidez 
lateral en unidades de T*m. La estructura sin refuerzo es la denominada (ESTi), 
con V invertidas (EST 2 ) y con X de multiples pisos (EST 3 ): 


KLEST 1 = 


27390.7 

-21825.3 

5152.3 

-21825.3 

37610.0 

-23723.3 

5152.3 

-23723.3 

37975.1 

-698.9 

5407.7 

-23765.2 

81.8 

-738.8 

5350.4 

51.5 

120.0 

-558.6 


-698.9 

81.8 

51.5 

5407.7 

-738.8 

120.0 

-23765.2 

5350.4 

-558.6 

37933.7 -23423.3 

4484.4 

-23423.3 35422.5 

-16679.5 

4484.4 

-16679.5 

12563.9 


KL_EST 2 = 


52583.5 

-39240.2 

5417.3 

-39240.2 

71910.9 

-41317.0 

5417.3 

-41317.0 

72221.2 

-534.2 

5538.1 

-41407.6 

169.0 

-664.4 

5441.4 

688.0 

631.0 

77.2 


-534.2 

169.0 

688.0 

5538.1 

-664.4 

631.0 

-41407.6 

5441.4 

77.2 

72125.7 

-41111.6 

5330.1 

-41111.6 

69585.7 

-33405.8 

5330.1 

-33405.8 

26658.0 
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52164.1 

-39659.2 

5956.3 

-707.7 

328.8 

1008.3 

-39659.2 

73265.8 

-41548.5 

5407.8 

-742.1 

126.5 

5956.3 

-41548.5 

71282.3 

-41593.6 

5681.6 

656.4 

-707.7 

5407.8 

-41593.6 

73589.5 

-41254.9 

4497.2 

328.8 

-742.1 

5681.6 

-41254.9 

68343.5 

-32342.2 

1008.3 

126.5 

656.4 

4497.2 

-32342.2 

26032.1 


Tambien se presenta la matriz de masas de las estructuras; para este 
caso se considera que es igual para los tres modelos. 



7.3128 

0 

0 

0 

0 

0 


0 

7.3128 

0 

0 

0 

0 

MASA = 

0 

0 

7.3128 

0 

0 

0 


0 

0 

0 

7.3128 

0 

0 


0 

0 

0 

0 

7.3128 

0 


_ 0 

0 

0 

0 

0 

4.142 _ 


En este punto se hace un parentesis al programa para generar el 
Pushover y con las matrices de rigidez y masa de cada estructura se puede 
generar un analisis modal espectral considerando el espectro de la NEC 15. El 
factor de reduccion de fuerzas sfsmicas R, se asume como 6; el tipo de suelo es 
D; la zona sfsmica V y la region Sierra, en base a ello se obtiene las propiedades 
dinamicas de cada estructura y se cumple con el proposito de calcular las derivas 
de piso, en la figura 40 se presentan mencionadas derivas. 


DERIVA POR PISO 



Figura 40 Derivas de piso de la respuesta modal espectral. 
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La deriva de piso maxima que tiene la estructura sin refuerzo es de 1.1% 
en los pisos inferiores, sin embargo, al emplear las diagonales de acero como 
mecanismos de reforzamiento esta deriva disminuye hasta 0.45%, 
aproximadamente, es decir que disminuye a mas de la mitad. No obstante, es 
importante notar que las diagonales en forma de V invertida o X de multiples pisos 
aportan una rigidez lateral similar a la estructura y por ello sus derivas son 
similares en cada piso, al igual que los valores en la matriz de rigidez. 

Los modelos de las diagonales de acero en forma de V invertida presentan 
la ventaja de generar espacios utiles en todos los pisos, mientras que las X de 
multiples pisos limitan la funcionalidad de espacio en los pisos que se colocan 
diagonales en forma de V. A pesar de esta ventaja, el emplear diagonales de 
acero en forma de V invertida en todos los pisos genera otros problemas 
relacionados con un desbalance de fuerzas en el nudo medio de la viga, formando 
un cortante vertical muy grande que debera ser absorbido por la viga y esta 
terminara fallando por flexion. Para contrarrestar el problema mencionado del 
desbalance de fuerzas se puede optar por diagonales en forma de X. 

A continuacion, en la figura 41 se muestra un fragmento del codigo para 
generar un modelo como el que se presenta en la figura 14, donde se divide a las 
vigas y columnas en 5 elementos, considerando los brazos rigidos, las zonas de 
plasticidad concentrada y los elementos elasticos. 


Codigo del programa principal, parte 6b: 


%% Elementos y nudos 

[Tabla_Elem_Nud,nod,X,Y]=ElementosNudos(nudt,nudcol,nod,NI,NJ,X,Y,d_col,d. 

beam); 

%. 

NI=Tabla_Elem_Nud(:,2)'; 

NJ=Tabla_Elem_Nud(:,3)'; 

[L,seno,coseno]=longitud(X,Y,NI,NJ); %long. de @elem, seno y coseno 
num_elem=nudt*5; 

HBuilding=max(Y); 

L_rig=0.6; 

NI2=[zeros(1 ,length(NI)), NI2(1 ,end-num_Diag+1 :end)]; 

NJ2=[zeros(1 ,length(NI)), NJ2(1 ,end-num_Diag+1 :end)]; 
[NI_Diag,NJ_Diag]=gn_portico3(NI,NJ, NI2,NJ2); 

%. 


Figura 41 Codigo del programa principal, parte 6b. 

8. MODELO CON ZONA DE PLASTICIDAD CONCENTRADA 
DE OPENSEES PARA DIAGONALES DE ACERO 

Este modelo considera plasticidad concentrada en las diagonales para el 
analisis Pushover no lineal en OpenSees, de manera similar al explicado en la 
seccion 4 del presente artfculo. Cada diagonal se divide en 5 partes como se 
muestra en la figura 42b, donde a los extremos se tiene elementos rigidos, luego 
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una rotula plastica de longitud nula, seguido de un elemento elastico de tipo 
armadura, es decir que trabaja a compresion y traccion. 



(b) 


Figura 42 a) Modelo generado en CEINCI-LAB b) Modelo para el analisis 
Pushover no lineal en OpenSees. 


En el modelo mostrado anteriormente se tiene 3 materiales en los 
elementos diagonales de acero; en la figura 43 se puede visualizar la ubicacion de 
las zonas de plasticidad concentrada. Estas zonas de plasticidad concentrada, de 
los elementos diagonales de acero, se formaran a una distancia L_rig. 

dool 


Rigid 

beam coiumn > 

• T. _ 

• > 

J V 


f— 

beam coiumn 


(b) 


♦ • 
4 


(a) 



Figura 43 a) Zonas de plasticidad concentrada en elementos diagonales de acero 
en la conexion viga - columna. b) Zonas de plasticidad concentrada en 
elementos diagonales de acero en el nudo medio de la viga. 
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Para el modelo de OpenSees se requiere nuevos nudos, coordenadas X, 
Y e identificacion de los nudos iniciales y finales que componen a cada elemento, 
a mas de ello es necesario identificar las secciones que se deben asignar a cada 
elemento. Se genera la funcion [Tabla_Elem_Nud_Diag, nod_Diag, X_Diag, 
YDiag] = ElementosNudosDiagonales(nod, numCOLVG, numDiag, 
NI_Diag, NJ_Diag, X, Y, L_rig); en este caso la matriz Tabla_Elem_Nud_Diag 
contiene la informacion relevante de nudos y elementos. Generando un modelo 
como se indica en la figura 44, donde se enumeran los elementos en las 
diagonales a partir de la ultima viga, es decir del elemento 300. 
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(a) (b) 

Figura 44 Numeracion de elementos; a) Porticos de diagonales en forma de V 
invertida; b) Porticos de diagonales en forma de X de multiples pisos. 

Este nuevo modelo se puede visualizar 60 elementos diagonales, es por 
ello que la matriz denominada Tabla_Elem_Nud_Diag tiene 60 filas donde cada 
una de ellas corresponde a un elemento. Esta matriz tiene 4 columnas donde la 
primera corresponde a una etiqueta, es decir numero de elemento original (figura 
36c y 36e) seguido de 2 ceros y el ultimo numero corresponde a los segmentos del 
1 al 5 en que fue dividido el elemento; la columna 2 y 3 corresponde al nudo inicial 
y final de cada nuevo elemento, la columna 4 es una etiqueta para identificar que 
tipo de propiedad corresponde a este elemento. El numero de la columna 4 puede 
ser 0 para elementos rfgidos (misma propiedad que un elemento elastico, pero con 
area e inercia multiplicada por 10), luego 1 para elementos inelasticos con 
comportamiento de armadura ( Truss ), 4 para zonas de plasticidad concentrada 
con propiedades de un elemento bilineal steel02. 
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Figura 45 Fragmento de Tabla_Elem_Nud_Diag a) Primera diagonal (elemento 61 
de la figura 36) b) Segunda diagonal (elemento 62 de la figura 36). 
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8.1. Modelo para zonas de plasticidad concentrada en las 
diagonales de acero 

El modelo para el comportamiento no lineal de las diagonales de acero 
denominado steel02 se basa en el modelo de Giuffre-Menegotto-Pinto. Este 
material es usado para la modelacion cfclica de acero estructural y varillas de 
refuerzo en hormigon armado. Sin embargo, debido a que el analisis es un 
Pushover monotonico, las diagonales de acero son simulados con modelos 
simplificados sin tener en cuenta el comportamiento real de estos elementos ante 
cargas reversibles. 

En la Figura 46 se muestra el comportamiento del material steel02. La 
curva constitutiva del material indica un comportamiento que se basa en una 
rigidez inicial, seguido de una zona de transicion exponencial hasta un 
endurecimiento del material (Bosco, Ferrara, Ghersi, Marino, & Rossi, 2014). 


▲ 



EP 


► 


£S = 0.002 

DeformaciCn unitaria (in/in) 

Figura 46 Modelo para las articulaciones en las diagonales de acero. 


9. TRANSFORMACIONES GEOMETRICAS EMPLEADAS EN 
OPENSEES 


Durante la definicion de los elementos vigas, columnas y diagonales de 
acero es necesario definir el tipo de transformacion de los parametros 
estructurales de rigidez y fuerzas resultantes. Las transformaciones geometricas 
se empelan para pasar del sistema de coordenadas fundamentales a un unico 
sistema de coordenadas globales como se puede visualizar en la figura 47. 

Elemento en sistema de Elemento en sistema de 

coordenadas globales . . coordenadas fundamentales 


V 


Fxi 



Mj 


► 


P 


Fyi 


Figura 47 Cambio de sistemas de coordenadas mediante transformacion 
geometrica. 
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Para explicar las transformaciones geometricas se emplea el vector de 
fuerzas de un sistema global denominado Q g i 0 bah relacionado con el vector de 
fuerzas de un sistema fundamental Qf un damentai■ 

Qfundamental \P ^j\ 06 ) 


Qglobal ~ [Fxi Fyi ^xj Fyj 


(17) 


La matriz de transformacion T, que relaciona los vectores de fuerzas como 
se indica en la ecuacion (18), se puede separar en dos transformaciones, una que 
realiza la transformacion de coordenadas fundamentales a locales y otra que 
realiza la transformacion de coordenadas locales a globales. 


Qglobal T ' Qfundamental 


(18) 


T — T a local ■ T a giobal (19) 

En el trabajo de Denavit y Hajjar (2013), se estudia el tema de las 
transformaciones geometricas, en el mismo se presenta las matrices de 
transformacion corrotacional al sistema de coordenadas global y local (matrices 20 
y 21), donde los parametros /? y L corresponden al angulo y la longitud de la barra 
como se muestra en la figura 48, los cuales con la transformacion corrotacional 
varfan en cada paso del analisis mientras que al emplear transformacion 
geometrica lineal estos valores son los iniciales /?o y Lo (matrices 22 y 23). 


_i 

l 


Figura 48 Geometria de un elemento tipo barra. 
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( 20 ) 


( 21 ) 


(22) 
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1 a local 


-1001 0 0 
0 1/Lo 1 0 —1/Lo 0 

L 0 1/Lo 0 0 1/Lo 1 


(23) 


Las matrices de transformacion P-Delta presentan la forma expuesta en 25 
y 26, donde el termino 8A y es la diferencia entre el desplazamiento transversal de 
los nudos, en el sistema local de coordenadas como se indica en la ecuacion 24. 
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0 

0 
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-1 

8A y /Lo 

0 
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0' 



T — 

1 a local 

0 

1/Lo 

1 

0 —1 /Lo 

0 





0 

1/Lo 

0 

0 1/Lo 

1 




(24) 


(25) 


(26) 


Matematicamente la matriz de transformacion lineal es constante, por lo 
tanto, se calcula unicamente al inicio del analisis con base a la geometrfa original 
de los elementos, es decir sin deformar. El calculo de la matriz de transformacion 
P-Delta requiere de un esfuerzo computacional adicional ya que tambien depende 
de los desplazamientos de los nudos del elemento, en cada paso del analisis, y 
mas aun la matriz de transformacion corrotacional, que debe ser calculada en 
cada paso del analisis en base a la geometrfa del elemento en dicho instante. En 
la tesis de Del Rincon (2017), se menciona que “la transformacion corrotacional 
toma en cuenta los cambios de longitud y posicion del elemento en cada paso 
mientras que la transformacion p-delta asume la teorfa de los pequenos 
desplazamientos”. 

OpenSees cuenta con estos tres metodos para realizar la transformacion 
geometrica de los elementos, es decir: lineal, P-Delta y corrotacional. 

En los modelos de las estructuras presentadas en este artfculo, debido a 
que el analisis es un Pushover monotonico y para optimizar el esfuerzo 
computacional se emplea la transformacion geometrica lineal para vigas, 
asumiendo piso rfgido (mismo desplazamiento horizontal en los nudos de las 
vigas); P-Detal para columnas porque basicamente se tiene desplazamientos 
laterales y la transformacion geometrfa corrotacional, en los elementos diagonales 
de acero, ya que las deformaciones son mayores y el cambio de su longitud es 
importante de considerar. 

10. DATOS DE OPENSEES GENERADOS POR CEINCI-LAB 


A continuacion en la figura 49, se muestra una parte del codigo en el cual 
se puede notar la funcion denominada Pushover_OpenSees_Diag(), que se 
presenta de color rojo para distinguirla, en este caso no existen datos de ingreso y 
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salida debido a que esta rutina trabaja con variables globales definidas al inicio del 
programa. Sin embargo, al ejecutarse esta linea, se genera un archivo con 
extension TCL cuyo nombre es Estructura_Pushover_Ref.tcl, este proceso es 
similar al que se visualiza en la figura 26. 


Codigo del programa principal, parte 7b: 


[Tabla_Elem_Nud_Diag,nod_Diag,X_Diag,Y_Diag]=ElementosNudosDiagonales(n 

od,num_COL_VG,num_Diag,NI_Diag,NJ_Diag,X,Y,L_rig); 

N l_D i ag=T abla_E le m_N ud_D i ag (:, 2)'; 

NJ_Diag=Tabla_Elem_Nud_Diag(:,3)'; 

[L_Diag,seno_Diag,coseno_Diag]=longitud(X_Diag,Y_Diag,NI_Diag,NJ_Diag); 
%long. de @elem, seno y coseno 
NI_Diag=[NI,NI_Diag]; NJ_Diag=[NJ,NJ_Diag]; 

%% Nudos con grados de libertad iguales 
[DOF]=MismosGDL(nod,np,Xcol,Ypisos,X,Y); 

Pushover_Opensees_Diag() 


Figura 49 Codigo del programa principal, parte 7b. 

Para verificar graficamente el modelo generado en CEINCI-LAB se puede 
elaborar los dibujos del portico en analisis, en este caso la funcion 
dibujoNudElemMODPUSH(X_Diag,Y_Diag,NI_Diag,NJ_Diag) presenta un 
grafico como se muestra en la figura 50b y 50c, en el cual se incluyen los nudos y 
elementos columnas, vigas y diagonales que componen al portico, con una 
distribucion similar a la presentada en la figura 44. 


Codigo del programa principal, parte 8b: 


%% Dibujos 

dibujoNudElemMODPUSH(X_Diag,Y_Diag,NI_Diag,NJ_Diag) 


Cvqwtms Cvtrwcturaf 




Figura 50 a) Codigo del programa principal, parte 8b. b) Estructura con diagonales 
en V invertida; c) Estructura con diagonales en X de multiples pisos. 
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El archivo Estructura_Pushover_Diag.tcl que se genera en la carpeta 
que se esta almacenando el programa en MATLAB, es ejecutado posteriormente 
por OpenSees y genera los calculos necesarios para el Pushover, la linea que 
indica esta accion se muestra a continuacion en la figura 51, donde el comando 
empieza con un signo de admiracion (!), luego se escribe OpenSees.exe seguido 
por el nombre del archivo tcl que describe el modelo, en este caso 
Pushover_Estructura_Ref.tcl 


Codigo del programa principal, parte 9b: 


! OpenSees.exe Pushover_Estructura_Ref.tcl 


Figura 51 Codigo del programa principal, parte 9b. 

Al ejecutarse OpenSees desde Matlab se generan los calculos del 
Pushover del modelo propuesto y aparece una ventana emergente con la 
simulacion del analisis como se muestra en la figura 52. 


(a) 
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Figura 52 Ventana con la simulacion del analisis Pushover en OpenSees. a) 
Estructura con diagonales en V invertida; b) Estructura con diagonales 
en X de multiples pisos. 

La funcion Graficar_Pushover2.m es la que genera la grafica con 
desplazamiento en el tope y cortante basal, es decir la curva del analisis Pushover, 
como se muestra en la figura 53b y 53c. En las graficas se indican las Ifneas que 
corresponde a los desplazamientos de fluencia y desplazamiento ultimo similar a 
la figura 31. 


Codigo del programa principal, parte IQb: 


%% Grafica Pushover 

Tfundamental=textread('Resultados_Modal\Periodo_Fundamental.out'); 
[Vbase,DT]=Graficar_Pushover2(nr,WBuilding,HBuilding,FP,Tfundamental); 
Resu ltados=[Vbase, DT]; 


(a) 
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Curva Pushover Curva Pushovor 




(b) 


(c) 


Figura 53 a) Codigo del programa principal, parte lOb. Curva Pushover generada 
por CEINCI-LAB mediante el acoplamiento con OpenSees. b) 
Estructura con diagonales en V invertida; c) Estructura con diagonales 
en X de multiples pisos. 


En la figura 54 se puede observar que las estructuras con diagonales 
aumentan la rigidez y resistencia ante cargas laterales. Las estructuras con 
diagonales presentan una ventaja, se puede apreciar un aumento de 
aproximadamente 40% en la fuerza maxima de corte en la base de la estructura. 


Curva Pushover 



0 0.25 0.5 0.75 1 

Desplazarmento en el tope. DT (m) 


Figura 54 Curva Pushover generada por CEINCI-LAB mediante el acoplamiento 
con OpenSees. 
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En el anexo 1 se presenta el codigo completo para el analisis no lineal 
estatico mediante la tecnica de Pushover, al inicio del programa se debe decidir 
para que estructura es el analisis. En este caso se entiende la variable 
ESTRUCTURA 1 para el portico sin reforzamiento, 2 para agregar los elementos 
en V invertida y 3 para diagonales en X de multiples pisos. 

11. COMENTARIOS Y CONCLUSIONES 

El objetivo de este artfculo es dar a conocer al lector las nuevas funciones 
de CEINCI-LAB y su acoplamiento a OpenSees, de modo que se facilite la 
modelacion y el ingreso de la estructura para su debido analisis. Es conocido que 
tanto CEINCI-LAB como OpenSees, pueden realizar independientemente el 
analisis Pushover, sin embargo, se busca ilustrar un enfoque aprovechando el 
potencial de ambos programas computacionales. 

La versatilidad del ingreso de los datos para la geometrfa de la estructura 
en CEINCI-LAB, resulta muy conveniente para el entendimiento del usuario y 
ahorra tiempo notablemente, el posterior analisis en OpenSees y la devolucion de 
resultados para procesar de la manera que considere mas adecuada el usuario, 
hacen de estas herramientas un conjunto provechoso para el investigador. 

Con el programa presentado en este artfculo, es posible generar un 
codigo para analisis Pushover de estructuras metalicas en OpenSees que puede 
ser modificado o se puede cambiar de forma parcial la programacion en funcion 
del criterio tecnico del usuario, lo que potencia un codigo flexible para nuevas 
investigaciones. 

Se muestra ademas el modelo de plasticidad concentrada, bajo el que se 
prepara el analisis y la forma de detallar cada elemento, para asl obtener los 
resultados de las curvas Pushover de cada estructura, es importante notar que se 
logra estimar una sustancial diferencia entre las propuestas de estructuras en 
analisis. 


Los mecanismos de reforzamiento propuestos son una buena alternativa 
cuando se requiere aumentar la rigidez lateral y resistencia de las estructuras 
ante solicitaciones de cargas laterales, como puede ser el sismo. Las diagonales 
en forma de V invertida presentan ventajas de arquitectura y funcionalidad, no 
obstante, se requiere evaluar que el desbalance de fuerzas en la viga no 
provoque una falla a flexion. Una alternativa para controlar este problema es 
colocar las diagonales en forma de X de piso a piso, esto a mas de aumentar la 
rigidez y resistencia elimina el problema del caso anterior. 
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ANEXO 1.0 


clear;clc 
format short 

% Acoplar OPENSEES con CEINCI-LAB 

global DOF nod sv sp nv np nr X Y Nl NJ n masa_nudo Xcol Ypisos L 
global coseno Es Abeam Ibeam Mybeam Ks_beamxL qp 
global Acol lcol Mycol Ks_colxL num elem Tabla_Elem_Nud 
global DP DN th_pP th pN th pcP th_pcN Dis_Fuerzas_Laterales T W fi 

%. 

global nod_Diag X_Diag Y_Diag NI_Diag NJ_Diag Abrace Ibrace num_Diag 

Tabla_Elem_Nud_Diag 

global Fy_brace Es_brace 

%% Geometrfa del portico 

Uso_Conexion = 1; % Si la conexion es Precalificada 1, sino 2 entonces Cpr = 1. 
ESTRUCTURA = 3; % 1 Sin refuerzo, 2 Con V invertida y 3 con X multiples pisos 
Graficar=0; % 0 si no quiere graficar, 1 en funcion del num de ESTRUCTURA, 
2 Todas 

sv =[7.315;7.315;7.315]; %lngresar un vector con la longitud de vanos en (metros) 
sp =[5.486;3.81; 3.81; 3.81; 3.81; 3.81]; %lngresar un vector con la altura de pisos 
en (metros) 

g=9.81; % Gravedad m/seg2 

qp=[-3.269;-3.269;-3.269;-3.269;-3.269;-1.8517]; % carga distribuida en vigas en 
cada piso T/m 

%. 

mar = [2;5;8;11 ;14;17]; % Para las V invertidas 
mar2 = [2;8;14]; % Para las V invertidas en el modelo con X 

marv = [5;11; 17]; % Para las V en el modelo con X 

%Para diagonales se calcula el area e inercia 
diae=0.18; %Diametro exterior 
diai=0.18-8*2/1000; %Diametro interior 
%. 

[nv,np,nudt,nudcol,nudvg,nod,nr]=geometria_nudo_viga(sv,sp); 

[X,Y]=glinea_portico2(nv,np,sv,sp,nod,nr); 

[NI,NJ]=gn_portico2(nr, nv, nudt, nudcol, nudvg); 

[CG,ngl]=cg_sismo2(nod,nr,Y); %Calcula los grados de libertad 
[VC]=vc(NI,NJ,CG); 

[L,seno,coseno]=longitud(X,Y,NI,NJ); %long. de @elem, seno y coseno 

num_COL_VG=length(L); 

Xcol=X(1,1 :nr); 

Ypisos=nonzeros(unique(Y)); 

%. 

[GEN]=geometria_nudo_diagonales(nv, np, nudt, mar); 
[NI2,NJ2]=gn_portico(GEN); 

[NI_Diag2,NJ_Diag2]=gn_portico3(NI,NJ, NI2,NJ2); 

%. 

[GEN]=geometria2_nudo_diagonales(nv, np, nudt,mar2,marv); 
[NI3,NJ3]=gn_portico(GEN); 

[NI_Diag3,NJ_Diag3]=gn_portico3(NI,NJ, NI3.NJ3); 

%. 
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NJ_Dib{1 }=NJ; NJ_Dib{2}=NJ_Diag2; NJ_Dib{3}=NJ_Diag3; 

%. 

if ESTRUCTURA ==1 
NI_Diag=NI; 

NJ_Diag=NJ; 

elseif ESTRUCTURA ==2 

NI_Diag=NI_Diag2; 

NJ_Diag=NJ_Diag2; 
elseif ESTRUCTURA ==3 
NI_Diag=NI_Diag3; 

NJ_Diag=NJ_Diag3; 

NI2=NI3; 

NJ2=NJ3; 

end 

[VC_Diag]=vc(NI_Diag,NJ_Diag,CG); 

[L_Diag,seno_Diag,coseno_Diag]=longitud(X,Y,NI_Diag,NJ_Diag); %long. de 
@elem, seno y coseno 
%% Masas 
W=zeros(np,1); 
for t=1 :np 
sum=0; 
for r=1 :nv 

sum=-qp(t,1)*sv(r,1); 

W(t,1)=W(t,1)+sum; 

end 

end 

sum=0; 

WBuiiding=0; 
for i=1 :np 

WBuilding=WBuilding+W(i,1); %Peso total en T 
for ij=1 :nv+1 
if ij==1 

NodalMass =(-qp(i,1 )/g)*sv(ij,1 )/2; 
elseif ij==nv+1 

NodalMass =(-qp(i, 1 )/g)*sv(nv, 1 )/2; 
else 

NodalMass =(-qp(i,1)/g)*sv(ij-1,1)/2+(-qp(i,1)/g)*sv(ij,1)/2; 
end 

masa_nudo(i,ij)=NodalMass; 

end 

MASA(i,i)=W(i,1)/g; 

end 

% Dis_Fuerzas_Laterales = 1; % 1 En funcion del primer modo de vibrar 
% 2 Distribucion parabolica; % 3 Distribucion triangular; % 4 Distribucion uniforme 
Dis_Fuerzas_Laterales = 1; 

% % Dibujo de la Estructura 2D 
if Graficar==0 
elseif Graficar==1 
dibujo(X,Y,NI_Diag,NJ_Diag) 
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dibujoNudElem(X,Y,NI_Diag,NJ_Diag,CG) 

dibujogdl(X,Y,NI_Diag,NJ_Diag,CG) 

else 

for op=1:3 

dibujo(X,Y,NI_Dib{op},NJ_Dib{op}) 

dibujoNudElem(X,Y,NI_Dib{op},NJ_Dib{op},CG) 

dibujogdl(X,Y,NI_Dib{op},NJ_Dib{op},CG) 

end 

end 

%% Propiedades del material 
Fy=25310.505; % Esfuerzo de fluencia, T/m2 
Es=200.0*1000000.0/9.81 ;% Modulo de elasticidad, T/m2 
Fy_brace=Fy; 

Es_brace=Es; 

Fu=40778.036; % Esfuerzo ultimo, T/m2 
if Uso_Conexion==1 

Cpr=(Fy+Fu)/(2*Fy); % Factor que toma en cuenta la resistencia maxima de la 
conexion. 

if Cpr < 1.2 

Cpr=Cpr; 

else 

Cpr =1.2; 
end 
else 
Cpr=1; 
end 

%% Secciones de los elementos 
% Column section % W30x173 
Zx = 0.0098; 

Acol = 0.0326; 
lcol = 0.0034; 

Mycol =Fy*Zx*Cpr; 
d_col= 773.2/1000.0; 

% Brace section % Tcircular Diametro de 18 cm y espesor de 8 mm 
Abrace=pi*diae A 2/4-pi*diai A 2/4; % area 
lbrace=pi*diae A 4/64-pi*diai A 4/64; % inercia 
% beam section % W30x99 
Zx = 0.0051; 

Abeam = 0.0186; 

Ibeam = 0.0017; 

Mybeam = Fy*Zx*Cpr; 
d_beam = 754.4/1000.0; 

% beam springs at Floor 2 # % W30x99 

DP= 0.886; % rate of cyclic deterioration for pos loading 

DN= 0.886; % rate of cyclic deterioration for neg loading 

th_pP= 0.0204; % plastic rot capacity for pos loading 
th_pN= 0.0204; % plastic rot capacity for neg loading 

th_pcP= 0.1001; % post-capping rot capacity for pos loading 
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th_pcN= 0.1001; % post-capping rot capacity for neg loading 

%. 

%% Propiedades dinamicas 
if ESTRUCTURA==1 


[Seccion]=Unica_Vg_COL(nudt,coseno,Acol,lcol,Abeam,lbeam); 
%Seccion=[ELEMENTO, BASE, ALTURA, NUMEROS DE ELEMENTOS A 
CREAR, INCREMETO EN NUMERO DE ELEMENTO] 
[ELEM]=gelem_portico(Seccion); 

% Matriz de Rigidez 
ELEMH=ELEM(1 :length(NI),:); 

[KH]=krigidez_acero(ngl,ELEMH,L,seno,coseno,VC,Es); 

KTT=KH; %<-MATRIZ "K" DE TODA LA ESTRUCTURA (Tn/m) 

else 

[Seccion]=Unica_Vg_COL(nudt,coseno,Acol,lcol,Abeam,lbeam); 
Seccion=[Seccion; length(NI)+1 Abrace Ibrace 11 1]; 
[ELEM]=gelem_portico(Seccion); 

% Matriz de Rigidez 
ELEMH=ELEM(1 :length(NI),:); 

[KH]=krigidez_acero(ngl,ELEMH,L,seno,coseno,VC,Es); 

%Contribucion de elementos Diagonales de Acero 
ELEMA=ELEM(length(NI)+1 :end,:); 

LA=L_Diag(length(NI)+1 :end); senoA=seno_Diag(length(NI)+1 :end); 
cosenoA=coseno_Diag(length(NI)+1 :end); VCA=VC_Diag(length(NI)+1 :end,:); 
num_Diag=length(LA); 


%%%%%%%%%%%%%%%%%%%%%%%%%%%% UNIDADES cm 


Est = Es/10; % Modulo de elasticidad, (kg/cm2) 

H=sp*100; 

hci=d_col*100;hcd=d_col*100;hdis=0;Lon=sv(2,1 )*100;hv=d_beam*100; 

Kef=zeros(np,1); 

for ie=1 :np 


[Kdiag(ie),Keq(ie)]=rigidez_equivalente_TADAS(hci,hcd,hdis,Lon,hv,H(ie),Kef(ie),A 

brace*10000,Est); 

end 

Kdiag=Kdiag/10; % Para tener en T/m 

ELEMT_Diag=[Kdiag(1 );Kdiag(1); 

Kdiag(2);Kdiag(2); 

Kdiag(3);Kdiag(3); 

Kdiag(4);Kdiag(4); 

Kdiag(5);Kdiag(5); 

Kdiag(6);Kdiag(6)]; 

[KA]=krigidez_tadas(num_Diag,ngl,ELEMT_Diag,senoA,cosenoA,VCA); 

%Rigidez de Diagonales 

KTT=KH+KA; %<-MATRIZ "K" DE TODA LA ESTRUCTURA (Tn/m) 

end 

na=np; 

kaa=KTT(1 :na,1 :na);kab=KTT(1 :na,na+1 :ngl); 
kba=kab'; kbb=KTT(na+1 :ngl,na+1 :ngl); 

KL=kaa-kab*inv(kbb)*kba; 

%. 

[T,fi,OM]=orden_eig(KL,MASA); 
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% Factores de participacion modal 

b=ones(np,1); 

for i=1 :np 

maxfi=max(abs(fi(:,i))); 
fi(:,i)=fi(:,i)/m axfi; 
end 

for i=1 :np 

FP(:,i)=abs(((fi(:,i))'*MASA*b)/((fi(:,i))'*MASA*(fi(:,i)))); 

end 

disp('Los periodos de vibracion en segundos con CEINCI-LAB son:') 

T % seg 

% Grabar matrices de masa y rigidez 

Masa_Rigidez = [MASA,KL]; 

mkdir Prop_Dinamicas_Estructuras_2D 

save(['Prop_Dinamicas_Estructuras_2D\EST_' num2str(ESTRUCTURA) 
'_M_K.txt'], ['Masa_Rigidez’],'-ascii'); %.DATO:Direccion carpeta 
% break 

%% Elementos y nudos 

[Tabla_Elem_Nud,nod,X,Y]=ElementosNudos(nudt,nudcol,nod,NI,NJ,X,Y,d_col,d_ 

beam); 

%. 

NI=Tabla_Elem_Nud(:,2)'; 

NJ=Tabla_Elem_Nud(:,3)'; 

[L,seno,coseno]=longitud(X,Y,NI,NJ); %long. de @elem, seno y coseno 
num_elem=nudt*5; 

HBuilding=max(Y); 

L_rig=0.6; 

NI2=[zeros(1 ,length(NI)), NI2(1 ,end-num_Diag+1 :end)]; 

NJ2=[zeros(1 ,length(NI)), NJ2(1 ,end-num_Diag+1 :end)]; 
[NI_Diag,NJ_Diag]=gn_portico3(NI,NJ, NI2,NJ2); 

%. 

[Tabla_Elem_Nud_Diag,nod_Diag,X_Diag,Y_Diag]=ElementosNudosDiagonales(n 

od,num_COL_VG,num_Diag,NI_Diag,NJ_Diag,X,Y,L_rig); 

N l_D i ag=T abla_E le m_N ud_D i ag (:, 2)'; 

NJ_Diag=Tabla_Elem_Nud_Diag(:,3)'; 

[L_Diag,seno_Diag,coseno_Diag]=longitud(X_Diag,Y_Diag,NI_Diag,NJ_Diag); 

%long. de @elem, seno y coseno 

NI_Diag=[NI,NI_Diag]; 

NJ_Diag=[NJ,NJ_Diag]; 

%% Nudos con grados de libertad iguales 
[DOF]=MismosGDL(nod,np,Xcol,Ypisos,X,Y); 

%% Dibujos 
if ESTRUCTURA==1 
dibujoNudElemMODPUSH(X,Y,NI,NJ) 

Pushover_Opensees() 

! OpenSees.exe Pushover Estructura.tcl 
else 

dibujoNudElemMODPUSH(X_Diag,Y_Diag,NI_Diag,NJ_Diag) 

Pushover_Opensees_Diag() 

! OpenSees.exe Pushover_Estructura_Ref.tcl 
end 
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%% Grafica Pushover 

Tfundamental=textread('Resultados_Modal\Periodo_Fundamental.out'); 

[Vbase,DT]=Graficar_Pushover3(nr,WBuilding,HBuilding,FP,Tfundamental); 
Resultados=[Vbase, DT]; 

%% Grabar los archivos 

mkdir Pushover_Estructuras_2D 

save(['Pushover_Estructuras_2D\EST_' num2str(ESTRUCTURA) 
'.txt'],['Resultados'],'-ascii'); %.DATO:Direccion carpeta 



